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Scalar Quarkonium Masses 

W.Lcc and D. Weingarten 

IBM Research, P.O. Box 218, Yorktown Heights, NY 10598 

We evaluate the valence approximation to the mass of scalar quarkonium for a range of different parameters. 
Our results strongly suggest that the infinite volume continuum limit of the mass of ss scalar quarkonium lies 
well below the mass of /j(1710). The resonance /o(1500) appears to the best candidate for ss scalar quarkonium. 
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For the valence approximation to the infinite 
volume continuum limit of the lightest scalar 
glueball mass, a calculation on GF11 Q] gives 
1740 ± 71 MeV, favoring /./(1710) as the lightest 
scalar glueball. Among other observed resonances 
which could have scalar glueball quantum num- 
bers, only / (1500) is near enough to 1740 ± 71 
MeV to be a possible alternate glueball candidate. 
Ref. Q suggests that /o(1500) is dominantly ss 
scalar quarkonium. Evidence for this identifica- 
tion is given in Ref. || . As a test of the interpreta- 
tion of /o(1500) as ss scalar quarkonium, we have 
now calculated, for Wilson quarks in the valence 
approximation, the mass of scalar quarkonium 
with two different values of lattice spacing and 
a range of different quark masses. Although the 
data we have is not sufficient to permit an extrap- 
olation of scalar quarkonium masses to the con- 
tinuum limit or to infinite volume, it is sufficient 
to show that the continuum infinite volume value 
for the valence approximation to the ss mass lies 
well below the mass of /j(1710). This result, 
combined with the discussion of Ref. ||, favors 
/o(1500) as ss scalar quarkonium and makes this 
interpretation for /j(1710) appear improbable. 

The calculations reported here were done on 
the GF11 parallel computer and required about 
5 months of operation at a sustained speed of ap- 
proximation 6 Gflops. 

To evaluate the mass of scalar quark-antiquark 
states, we define operators for these states by 
first fixing each gauge configuration to lattice 
Coulomb gauge, then convoluting, in the space 
direction, the Wilson quark field ^ sc (x) with a 
gaussian G r {x) with root-mean-squared radius r 
to produce a smeared Coulomb gauge quark field 



^ r sc ( x ) (§■ Here 1 < s < 4 and 1 < c < 3 are spin 
and color indices, respectively. We assume only a 
single flavor of quark. In a gamma-matrix repre- 
sentation with 75 given by the diagonal entries 1, 
1, -1, -1, define the upper and lower quark fields 
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and define W^(x) and ^> r c d (x), similarly from 
^ sc (x) with £* in place of £ s . Here the £ s , for 
each s, are independent random cube roots of 1. 
The scalar operator 

S£(x) = W u (x)^ ru {x)+W d {x)^ rd {x), 
when averaged over £ becomes the more familiar 
S r (x) = W(x)^ r (x). 

A pseudoscalar P£(x) and its average P r {x) can 
be defined in analogy to S£(x) and S r (x). 
For scalar correlation functions we choose 

C r s (t) = ^2<S r (x,t)Sl(0)>- 

X 

J2<s r (z,t)><si(o)> 

where averages are over £ and gauge field con- 
figurations, with one random vector (£i,...£4) 
for each field. A pseudoscalar C T P (t) can be 
defined similarly. The correlations C r s {t) and 
Cp(t) together require a factor of six fewer quark 
matrix inversions per gauge configuration than 
needed for C'g{t) and C£{t) defined with S r (x) 
and P r {x), respectively, used as both source and 
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sink operators. To obtain propagators with a 
fixed statistical uncertainty, however, some of 
this factor of six will be lost to the additional 
noise arising from £. Figure |l| shows the actual 
gain in arithmetic work for the scalar correlation, 
6[D' s (t)/D s (t)] 2 , where D s (t) and D' s (t) are the 
statistical dispersions in C r s (t) and C'g (t) respec- 
tively. Figure |^ shows the corresponding actual 
gain in arithmetic work for the pseudoscalar cor- 
relation. The values in Figures [l] and |^ were found 
from 188 independent gauge configurations on a 
lattice 16 3 x 24 with (3 of 5.70 and k of 0.1650. 
For smaller k we expect the gain to be greater. 
For the time intervals of our scalar mass fits, the 
actual gain turns out to be about a factor of 2. 
Without this factor, our five months of calcula- 
tion would have required ten. 

The cost of generating gauge field configura- 
tions and gauge fixing we found could be reduced, 
for the scalar correlation function, by evaluating 
propagators for several different source time val- 
ues on each lattice. 

Figure [| shows effective masses, the fitted mass, 
and the fitting range for the scalar correlation 
function given by 1972 gauge configurations on 
a lattice 16 3 x 24 with (3 of 5.70, k of 0.1650 and 
smearing radius r of y6- For each gauge configu- 
ration, propagators were found from six different 
starting times. Figures [| shows effective masses, 
the fitted mass and the fitting range for the scalar 
correlation function given by 1733 gauge configu- 
rations on a lattice 24 4 with (3 of 5.93, K of 0.1567 
and r of y/27/2. For each gauge configuration, 
propagators were found from four different start- 
ing times. In physical units, the quark mass and 
smearing radius in Figures || and |] are nearly 
equal. Scalar and pseudoscalar masses for several 
different k at each f3 are given in Tables [l] and 
||. In Ref. [Q, the rho mass in lattice units for 
16 3 x 24 with /3 of 5.70 and for 24 4 with (3 of 5.93 
was reported to be 0.5676(79) and 0.3851(79), re- 
spectively. Thus in physical units the lattice pe- 
riod at the two different values of j3 we consider 
are nearly equal. 

The ss scalar mass we found by interpolation 
in quark mass to the strange quark mass cor- 
responding to hopping constants determined in 
Ref. f§ to be 0.16404 at (3 of 5.70 and 0.15620 
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m sc 


range 


x 2 


0.1625 


1.299(12) 


3 - 5 


0.320 


0.16404 


1.291(10) 






0.1650 


1.287(6) 


2 - 4 


0.002 
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7Tip S 


range 


x 2 


0.1625 


0.5795(6) 


7 - 10 


0.371 


0.1650 


0.4560(5) 


7 - 10 


0.250 



Table 1 

Masses, fitting ranges and \ 2 P er degree of free- 
dom of mass fits on a 16 3 x 24 lattice at (3 of 
5.70. 



K 


m sc 


range 


x 2 


0.1539 


0.856(3) 


4 - 


11 


1.39 


0.1554 


0.806(4) 


4 - 


11 


1.40 


0.1562 


0.788(17) 








0.1567 


0.777(5) 


4 - 


10 


1.44 
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TYlps 


range 


x 2 


0.1539 


0.4819(5) 


8 - 


12 


0.62 


0.1554 


0.3982(5) 


8 - 


11 


0.26 


0.1567 


0.3147(6) 


8 - 


11 


0.52 



Table 2 

Masses, fitting ranges and \ 2 P er degree of free- 
dom of mass fits on on a 24 4 lattice at (3 of 5.93. 

at (3 of 5.93. The resulting masses, in units of 
trip 0] , are shown in Figure |as a function of lat- 
tice spacing in comparison to the predicted scalar 
glueball mass and the masses of /j(1710) and 
/o(1500). Figure || stongly suggests that for the 
fixed volume we are using, the continuum value 
of the ss scalar mass will lie well below the mass 
of /j(1710). Since the ss scalar has one unit 
of orbital angular momentum, its physics radius 
should be larger than that of the ss <f> meson with 
no orbital angular momentum. The scalar mass 
divided by should therefore fall with volume. 
Volume dependence data for in Ref. |Q] then 
implies that infinite volume values of the scalar 
mass divided by m p will be at most 2.8% above 
the numbers shown in Figure ^[ Our qualitative 
conclusion about the scalar ss mass remains the 
same in infinite volume. Thus /o(1500) appears 
to be a likely candidate for the the ss scalar, while 
this assignment for /j(1710) appears quite un- 
likely. 
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Figure 1. For the scalar propagator on a 16 3 x 24 
lattice with (3 of 5.70 and k of 0.1650, the perfor- 
mance gain of a random source in comparison to 
a causal source. 
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Figure 3. Effective scalar mass, fitted mass and 
fitting range for a 16 3 x 24 lattice with (3 of 5.70 
and k of 0.1650. 
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Figure 2. For the pseudoscalar propagator on a 
16 3 x 24 lattice with of 5.70 and n of 0.1650, 
the performance gain of a random source in com- 
parison to a causal source. 
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Figure 4. Effective scalar mass, fitted mass and 
fitting range for a 24 4 lattice with (3 of 5.93 and 
k of 0.1567. 
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Figure 5. Scalar ss quarkonium masses for two 
different values of lattice spacing. 



